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Abstract 



We present waveform observations of electromagnetic lower hybrid and whistler waves with 
f, <C f < f„ downstream of four supercritical interplanetary (IP) shocks using the Wind search coil 
magnetometer. The whistler waves were observed to have a weak positive correlation between 
SB and normalized heat flux magnitude and an inverse correlation with T t , ; ,/T„. All were observed 
simultaneous with electron distributions satisfying the whistler heat flux instability threshold and 
most with Tjj/Tp > 1.01. Thus, the whistler mode waves appear to be driven by a heat flux 
instability and cause perpendicular heating of the halo electrons. The lower hybrid waves show a 
much weaker correlation between <5B and normalized heat flux magnitude and are often observed 
near magnetic field gradients. A third type of event shows fluctuations consistent with a mixture 
of both lower hybrid and whistler mode waves. These results suggest that whistler waves may 
indeed be regulating the electron heat flux and the halo temperature anisotropy, which is important 
for theories and simulations of electron distribution evolution from the sun to the earth. 



1 Introduction 

Electromagnetic instabilities have been postulated to occur in and around collisionless shock waves 
for over 40 years [ Tidman and Krall 1971) . A list of the waves/instabilities that are most commonly 



proposed to provide significant anomalous resistivity can be found in Wu et al. [ 1984| . Two types 
of waves have been of particular interest in collisionless shock energy dissipation: lower hybrid and 
whistler waves. Lower hybrid waves (LHWs) is a term used to encompass a broad range of waves 
and/or instabilities in plasma physics. Though LHWs are often thought to exist as a simple linearly 
polarized electrostatic wave propagating perpendicular to the ambient magnetic field with f < f ;/ , (= 
Vfafce for solar wind parameters), they can exist as a right-hand circularly polarized electromag- 
netic mode | Wu et a/.] 1 1983) . Measurements at the terrestrial magnetopause have found that LHWs 



not only have kii ^ 0, the electromagnetic component can be comperable to the electrostatic one 
jMarsch and Chang\\m3\\Silin et a^||2005| . 
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Whistler mode waves (WWs) can exist as a right-hand polarized electromagnetic mode prop- 
agating either parallel or at oblique angles to the magnetic field or as a quasi-electrostatic mode 
propagating at highly oblique angles with frequencies f cj < f < f re |gnce|[T964[ . In the solar wind 
at low frequencies (f r , <C f < f//,), whistler mode waves can be very oblique electromagnetic modes 



which can lie on the same branch of the dispersion relation as LHWs \Marsch and Chang 1983 1 



and/or magnetosonic modes | Wu et al. 1983]. Higher frequency (f lh < f <C f ce ) WWs for which f — > 
f ce have k ± — ► \ Kennel and Petscheck 1966] and can couple to multiple wave modes \Dyrud and 
\Oppenhelm\ |2006 1 . 

Both LHWs and WWs can result from a number of different free energy sources. WWs can 



be unstable to the heat flux carrying electrons [e.g. Gary et al. 1994| [19991 and/or temperature 



anisotropy [e.g. \Kennel and Petscheck 1966) in the solar wind. The relevant LHW free energy 
sources at collisionless shocks include, but are not limited to, currents [e.g. Lemons and Gary||1978| 
\Matsukiyo and Scholer 2006|, ion veloc ity rings [e.g. \Akimoto et a/.|[T985) , and/or heat flux car- 



rying electrons [e.g. Marsch and Chang 1983). More importantly, both modes can interact with 
electrons and ions and exchange energy between the two species. It is well known that LHWs 
and highly oblique WWs can interact with high energy electrons parallel to the magnetic field and 
the bulk of the ion distributions perpendicular to the magnetic field through Landau resonance and 
nonlinear trapping | Savoini and Lembege 1995 Wu et al. 1983| . This is due to their large phase 
velocity parallel to the magnetic field (co/ky) and much smaller perpendicular phase velocity (co/k ± ). 
The stochastic energization of electrons and ions results in a nonthermal tail of the electrons parallel 
to the magnetic field and perpendicular ion heating. Observation of these features in the ion and 
electron distributions in the ramp of one supercritical IP shock has been attributed to highly oblique 
whistler mode waves | Wil son III et al. \ [20 1 2 1 . 



Although there have been numerous studies of LHWs in the magnetosphere [e.g. Cattell et al. 



1995 1, there have only been a few observations done at the terrestrial bow shock | Mellott and Green- 



stadt 



1988 



Walker et al. 2008 Wygantetal. 1987 1 and only one observational study at an IP shock 



| Zhang and Matsumoto 1998 ] . WWs have been observed in nearly every space plasma environment 
including, but not limited to, the terrestrial magnetosphere [e.g. Russell et al.\ |1969), in the solar 
wind [e.g. Coroniti et al. 1982) , upstream of IP shocks [e.g. 



Russell et al. 1983 Wilson III et al. 



2009], downstream of IP shocks \Moullard et al. 1998 2001 1, upstream of planetary bow shocks 
[e.g. Hoppe et al. 1981], and in the ramp of supercritical IP shocks | Wilson III et al. 2012 1. 



Electron distributions in the solar wind are composed of a cold dense core (subscript c), a hotter 
more tenuous halo (subscript h), and a narrow field-aligned strahl [Feldman et al. |1973). Many the- 



ories [e.g. \Gary et al\ 1 1 994) \Saito and Gary\ |2007| |Vbcfa et al\ |2005) of the evolution of electrons 
from the Sun to the Earth suggest WWs scatte r strahl electrons into the halo, which may explain 
the observed changes [e.g. Stverdk et al. 2009) in the relative densities of strahl versus halo elec- 
trons. Theories have also suggested that WWs constrain the heat flux carrying electrons [e.g. \Gary 



et al. 1994 1999] and the halo/core electron temperature anisotropies [e.g. Saito et al. 2008; Vocks 



and Mann 2003 



[e.g. Salem et al. 



The core electrons have been observed to isotropize due to binary collisions 



2003], while halo electrons appear to be constrained b y the whistler anisotropy 



and fireshose instabilities in the slow solar wind [e.g. Stverdk et al. [2008J . These studies did not, 
however, examine waveform captures at the time of the particle distributions. Thus, simultaneous 
observations of particle distributions with waveforms can test theories and improve our understand- 
ing of wave-particle dynamics in the solar wind. 

In this paper, we present an analysis of 47 individual waveform captures, sampled at 1875 sam- 
ples per second, observed by the Wind search coil magnetometer downstream of four supercritical 
IP shocks. The shocks were selected because they had search coil waveform captures with power 
peaked at f„ <C f JC < f/;,, where f sc is the spacecraft frame frequency. The particle distributions ob- 
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served simultaneously satisfy the whistler heat flux and whistler anisotropy instability thresholds of 
\Gary et aZ.|fT994) <md \Gary et g£] | fT999l 

The paper is organized as follows: Section [2] introduces and outlines the data sets and analysis 
techniques; Section [3] describes the waveform observations; Section |4] explains how we estimate 
wave numbers using multiple methods; Section [5] discusses our analysis of possible free energy 
sources and interpretation of wave modes; and Section [6]provides conclusions of our study. 



2 Data Sets and Analysis 



Waveform captures were obtained from the Wind/WAVES instrument \Bougeret et al. 1995 1, using 



the time domain sampler slow (TDSS) receiver, which provides a waveform capture (herein called 
TDSS event) of 2048 points with timespans of ^1000 ms in this study with a lowest resolvable 
frequency of ^3 Hz. Therefore, waves with spacecraft frame frequencies below 3 Hz will not be 
observed. For all the events in this study, TDSS events are composed of three magnetic (8Bj) and 
one electric {8Ej) field component, where j = x, y, or z. We use 5 to distinguish between the 
high frequency fields measured by the TDS instrument and the quasi-static magnetic fields mea- 
sured by the fluxgate magnetometer. Dynamic waveform analysis was performed through the use 
of the Morlet wavelet transform \Torrence and Compo\ |1998a[ for each component of the TDSS 
events examined in this study. For every TDSS wavelet presented herein, the spectra from all three 
components will have the same color scale range. The Wind/WAVES instrument also contains an 
onboard time-averaged spectral intensity instrument, the thermal noise receiver (TNR), used to de- 



termine the local electron density (n,,) from the plasma line | Meyer- Vernet and Perche 1989 1. For 
more details about the WAVES instrument see | Bougeret et al. 1995| , and for analysis see Wilson III 
et al. | 2010| . High time resolution magnetic field data was obtained from the dual triaxial fluxgate 



magnetometers \Lep ping et al.\\\995\ . The time resolution is ^0.092 s or ~11 samples/second. 



We obtained full An steradian low energy (<30 keV) electron and ion distributions from the 
Wind/3DP EES A and PES A particle detectors |Lw et a/.|[T995) . Both EES A Low and PES A High 
can return full An steradian distribution functions once every spin period (~3 s) in burst mode. 
There was burst mode particle data for the 1998-08-26 and 2000-02-1 1 events but not the 1998-09- 
24 event. The EES A Low detector measures electrons with < 1 . 1 keV. We use EES A Low (electrons) 
and PESA Low (ions) to calculate moments of the distribution functions and the anisotropy factor, 
A, = Txj/Tii j - 1, where j corresponds to the jth-component of the electron distribution (e = entire 
distribution, c = core, and h = halo) and _L(||) are defined with respect to ambient magnetic field, 
B„. The method for determining the break energy between halo and core electrons and the method 
to calculate the electron heat flux are outlined by Wilson III et al. [ 2009 1 . 

Supplemental ion distributions were calculated from the two Faraday Cup (FC) ion instruments 
from the Wind SWE experiment [Ogilvie et al. 1995}. The SWE FCs can produce reduced ion 



distribution functions with up to 20 angular and 30 energy per charge bins every 92 seconds over an 
energy range of ^150 eV to ^8 keV (^1200 km/s proton) | Kasper et al. 2006) . We will use the 
results of the reduced distributions produced by SWE to compare with our results from 3DP. 

The relevant shock parameters, determined by Kasper [2007], for the events studied are shown 
in Table[T] where n is the shock normal vector, V j; ,„ is the shock normal speed in the spacecraft (SC) 
frame, 6 B „ is the angle between the shock normal vector and the upstream magnetic field vector, My 
is the fast mode Mach number, and N i2 /N n is the shock compression ratio. The first critical Mach 



number, M„ , was estimated for each IP shock using the methods outlined by Edmiston and Kennel 



[ 1984| assuming a polytrope index, y = 5/3. All of the shocks examined in this study were found to 
have M f /Mcr > 1 (last column of Table[T]i. 
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The wave vector, k, and the polarization with respect to the ambient magnetic field were deter- 
mined using Minimum Variance Analysis (MVA) \Khrabrov and Sonnerup 1998 1. The low (f sc j on .) 
and high (f, c ,;„ g ;,) frequencies for each bandpass filter are determined from SC frame wavelet analysis. 
Multiple standard Fourier bandpass filters were applied to each waveform and then we used MVA 
on specific subintervals. The bandpass filters were applied to the entire TDSS event to reduce edge 
effects in the analysis of each subinterval. Wave vector solutions are kept if A,„ w /A m ,„ > 10, where 
Xj is the ] th eigenvalue of spectral matrix with X max > X mld > X mi „. The details of this technique are 
discussed m \Wilson III et al.\$2M)\ . 

The TDS data are displayed in two despun coordinate systems, GSE coordinates and field- 
aligned coordinates (FACs). We define the FAC system with one axis parallel to B fl , a second parallel 
to B /|B„| x n, and the third axis completes the right-handed system. The long duration of the TDSS 
events with respect to spacecraft and variability in the magnetic field direction resulted in our us- 
ing the instantaneous magnetic field interpolated from the fluxgate magnetometer measurements to 
rotate each vector in every TDSS event. 



3 Waveform Observations 

In Figure[T]we show examples of each of the three waveform types to be examined in this study and 
summarize their physical properties. We will explain how these properties were determined in more 
detail in the following section. The three FAC components (5Bii (blue); 8B ±A (orange); and 8B ±2 
(green)), and a corresponding Morlet wavelet transform are plotted. Overplotted on each wavelet 
transform are the lower hybrid (LH) resonance frequency, f ft , and the electron cyclotron frequency, 

fee- 

Figure [IJV is an example of an electromagnetic LHW. The properties used to identify these 
modes in this study are: (1) irregular waveform with power peaked at f JC < f lh ; (2) mixtures of right 
and left-handed elliptical polarizations (SC Frame) with respect to the magnetic field; (3) typically 
propagate at angles greater than 45° from the magnetic field; (4) propagation direction shows no 
apparent dependence on solar wind direction; (5) have no apparent fluctuations for f sc > 60 Hz; (6) 
exhibit significant amplitudes in all three FAC components; and (7) typically have amplitudes >2 
nT peak-to-peak. There were 23/47 TDSS events classified as electromagnetic LHWs downstream 
of the 1998-08-26, 1998-09-24, and 2000-04-06 shocks. 

The WWs (Figure[TJi) are characterized by six properties: (1) peak power at f sc > 60 Hz; (2) 
smaller amplitudes than the electromagnetic LHWs (~1 nT peak-to-peak); (3) near circular right- 
hand polarization (SC Frame) with respect to the magnetic field; (4) only small 8Bf (5) typically 
propagate within 30° of the magnetic field; and (6) typically propagate at >45° from the solar wind 
velocity. There were 13/47 total TDSS events classified as WWs downstream of the 1998-08-26 and 
2000-02-11 shocks. 

The third type of waveform observed in our study (Figure [T]C) is a mixture of both LHWs and 
WWs. These waves (MIXED) are characterized by five properties: (1) the majority have large am- 
plitudes (>2 nT peak-to-peak) for f sc < f„, and smaller amplitudes (~1 nT peak-to-peak) for f sc > f„,; 
(2) right-hand polarizations (SC Frame) in both high and low frequencies; (3) significant amplitudes 
in SB ii, 5B ±1 , and 8B ±2 for low frequencies but only 8B ±A and 8B ±2 for higher frequencies; (4) 
lower and higher frequency components show propagation characteristics consistent with the LHW 
and WWs, respectively; and (5) observed near sharp magnetic field gradients. We observed 1 1/47 
TDSS events classified as MIXED downstream of the 1998-09-24, 2000-02-11, and 2000-04-06 
shocks. 

At first glance, these modes appear to be the simultaneous occurrence of LHWs with WWs su- 
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perposed. However, we will show in Section 5. 1 that the MIXED modes are observed under plasma 
conditions that are different enough from those found with the LHW and WWs to warrant a distinc- 
tion. In Section |4~T| we will define the separation between the high and low frequency components 
for the MIXED modes and discuss this in more detail. 



4 Wave Number and Rest Frame Frequency Estimates 

In this section, we first discuss and summarize our minimum variance results in Section |4~T| Then 
we discuss and summarize our estimates of the wave numbers and rest frame frequencies in Section 

4.1 Minimum Variance Results 

As discussed in Section [2] we perform MVA on multiple frequency filters for multiple subintervals 
of each TDSS event. Each subinterval yielding X mid lX mi „ > 10 gives us an estimate of the wave 
vector (k), which we can use to calculate wave propagation angles between the ambient magnetic 
field (9 kB ) and local V. IK (9 kv ). Note that for single point measurements with only magnetic fields, 
the sign of the wave vector cannot be definitively determined. Therefore all angles are normalized 
to a 0° to 90° scale. 

Table [2] presents a summary of the wave observations where the first column defines the wave 
type, the second column shows the number of TDSS events for each wave type, the third column 
shows the number of wave vectors satisfying the conditions discussed above, and the fourth column 
shows the mean plus or minus the standard deviation of the mean wave amplitudes. Note that the 
wave amplitudes were determined from the filtered fields for each wave type. 

We examined 47 TDSS events and determined a total of 461 unique wave vectors from MVA. 
-70% of the 118 LHW wave vectors satisfied 9 kB > 45° and -87% had 9 kB > 30°. The LHW 
wave vectors showed a broad range of kv peaked near 45°. Nearly all of the 138 WW wave vectors 
satisfied 9 kB < 45° and -83% had kB < 30°. All of the WW wave vectors had 9 kv > 45° and 
-91% had G kv > 60°. 

For the MIXED waves, the high and low frequency components were separated before wave 
normal angles were determined. The TDSS events in our study all occurred under conditions with 
f/;, < 40 Hz. In the SC frame, the power is primarily at frequencies <40 Hz and/or >60 Hz. We 
define high frequency (f;„ ; ;,) in the SC frame as (f sc > f //( ) and (f sc > 40 Hz). We define low frequency 
(f w ) in the SC frame as (f IC < f„,) or (f sc < 40 Hz). We found 43/205 MIXED wave vectors satisfied 
f tow and had a broad range of 9 kB and 9 kv , both peaked near 45°. 162/205 MIXED wave vectors 
satisfied f Kgh and -67% had 9 kB < 45° and -78% had 9 kv > 45°. In summary, the MIXED modes 
satisfying f high have similar 9 kB and 9 kv as the WWs and those satisfying f lmv are similar to the LHWs. 
While their propagation characteristics are consistent with the LHWs and WWs, they are observed 
under different plasma conditions than the LHWs(WWs). 

4.2 Doppler Shift Results 

In the solar wind, waves can propagate at oblique angles with respect to the bulk flow (i.e. 9 kv ^ 
90) causing a relatively stationary detector to observe a Doppler-shifted frequency. In the spacecraft 
(SC) frame of reference, the wave frequency is a Doppler shifted signal given by Equation [3] To 
determine the rest frame frequency using single point measurements we need to assume a wave 
dispersion relation. In this section, we will compare our wave number estimates from four different 
methods for one example (described in Appendix|A|) and then discuss the results for all TDSS events 
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analyzed. 

The results obtained for the waveform in Figure [T]C calculation are shown in Table [5] We used 
two frequency filters on this wave (f sc ~ 5-30 Hz and ~ 120-200 Hz) and so Table[3]is separated into 
two sections. The first part of each section of the table shows the background plasma parameters 
for this event along with the filter frequency, 9 kB , and kv ranges. The second part of each section 
shows the wave number estimates. 

Table [3] shows an increase in kv with increasing f sc for the example wave. We found this to 
be a common trend with all the waves examined herein. Although we do not have a definitive 
explanation why 9 kv increases with increasing f w , it is possible that waves with small kv at higher 
frequencies are Doppler shifted above the Nyquist frequency of the search coil magnetometer (^940 
Hz). Note that the largest value of f sc ~ 400 Hz, which is less than half the Nyquist frequency. In 
our analysis, we use Equation]?] for both the WWs and MIXED modes satisfying f Mgh and Doppler 
shifted Equation[6]for LHWs and MIXED modes satisfying f law . 

A summary of the results are shown in Table|4]for kc/(O IK , kp„, k||p re , k ± p«, and co/Cl. Here, 
= (On, (= V^A) for the LHWs and the MIXED modes satisfying f lm and D. — i2„ for the 
WWs and MIXED modes satisfying f Ugh . For each wave type there are two rows, where the first row 
shows the absolute ranges and the second shows the mean plus or minus the standard deviation of 
the mean. 

The results for the LHWs and MIXED modes satisfying f !aw show a broad range of kc/(O pe and 
£0/0;/,. The broad range of f sc may be due to Doppler broadening resulting from a range of 9 kv for 
each TDSS event. These results are consistent with theory \ Marsch and Chang\ |1983[ \Wu et al. 
1983 1 . The broad range of (o/£l ce and kc/o^e for the WWs and MIXED modes satisfying f/ 7I ^ are 



consistent with theory [e.g. Kennel and Petscheck 1966 1 and previous observations \Coroniti et al. 



\\9%2\\Lengyel-Frey et al.\\\996\\Zhang and M atsumoto 1998 1 of whistler mode waves. Thus, we 
argue that the observed modes are indeed on the whistler branch of the dispersion relation. 

The frequency spectrum, oblique propagation, and amplitudes of the LHWs are consistent with 
previous observations at collisionless shocks \Walker et al.\ [20 08 ; Zhang and Matsumoto, 1998 1. 
The waveforms also show a resemblance to previous observations of a highly oblique class of waves, 
identified as whistler mode waves, observed downstream of the bow shock \Zha ng et al.\ 1998 1. 
However, the LHWs presented herein are much larger in amplitude and observed much farther from 
the shock ramp. The near circular right-hand polarization, amplitudes, and observed frequencies 
of the WWs are consistent with previous observations of narrow band whistler mode waves at IP 
shocks \Moullard et al.\ 1998 200T Zhang and Matsumoto 1998 1 and stream interaction regions 



\Breneman et al. 2010 



At higher frequencies the wave vectors become more aligned with the 
magnetic field and more oblique to the solar wind velocity, consistent with theory [e.g. Kennel\ 



and Petscheck 1966 1 and previous observations \Lengyel-Frey et al. \996\\Zhang a nd Matsumoto 
1998]. 



5 Free Energy and Interpretation of Waves 

In this section, we first introduce some previous work and theory on electromagnetic LHW and 
then WWs in the solar wind. Following this, we will discuss the observed electron temperature 
anisotropies and normalized heat flux magnitude. Finally we will discuss possible ion free energy 
sources. 
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5.1 Possible Free Energy Sources 

There are multiple theories for the free energy sources for waves with f„ < f < f ce in the solar 
wind. Gary et al. | 1994 1999); Marsch and Chang |1983| proposed that waves with f ci <C f < 
f ce can extract free energy for wave growth from the heat flux carrying electrons in the solar wind. 
Electromagnetic LHWs may also be driven by currents and/or gyrating ions \Akimoto et al. 1985 1. 



The theories for WW generation involve electron temperature anisotropies and the electron heat flux. 
Kennel and Petscheck | 1966 1 showed for a single bi-Maxwellian distribution that if: 

co ^ A e 



(1) 



was satisfied, the distribution would be unstable to the whistler mode. Gary et al. | 1994J showed 



that, even if no finite electron heat flux exists, if T«f/Y« c is small and T ±;i /T| A > 1.01, then the whistler 



anisotropy instability may be excited. Furthermore, Gary et aT\ \ 1999 1 determined that when T xh fTp 
> 1.01 in the presence of a finite electron heat flux, the heat flux carrying electrons are always un- 
stable to whistler heat flux instability. 

2009) , we found 



In an earlier study of whistler mode waves near IP shocks | Wilson III et al. 



that the waves were observed simultaneously with particle distributions unstable to the whistler heat 
flux instability. Based on estimates of the parameters of Gary et al. 1 1994) , all of the WWs satisfied 
the threshold criteria for a whistler heat flux instability. Furthermore, ~77% (10/13) of the WWs 
were observed when T ±h /T^ h > 1.01 and with a finite electron heat flux. Interestingly, 19/23 of the 
electromagnetic LHWs were observed with T ±h /T« h > 1.01 and satisfied the same threshold criteria 



of Gary et al. | 1994 1 as the WWs. For the MIXED modes, however, only 6/1 1 satisfied the threshold 
criteria of Gary et al. [ 1994|. 

Figure [2] shows the filtered peak-to-peak wave amplitudes, 8B, plotted vs. T eh /T ec (panels A and 
C) and normalized heat flux magnitude (panels B and D), |q e |/q , for the WW and LHWs. Recall 
that for each TDSS event, we used multiple frequency filters and examined multiple subintervals, 
each with their own peak-to-peak 8B. Thus, the red squares indicate the mean values of 8B and the 
error bars show the range of values for each TDSS event for all filters and subintervals. We used 
green dashed lines to highlight qualitative trends in the data. Note that the MIXED modes showed 
no relationship between SB and any of the parameters examined in our study (not shown). 

Figures |2ji and[2)) show a weak positive correlation between 8B and |q £ |/q„. One can see that 
the correlation is much stronger for the WWs than the LHWs and that the relationship for the LHWs 
only exists for |q f |/q D > 0.03. Figure|2jV shows a weak inverse relationship between SB and T eh /T ec 
for the WWs, which may be explained by increased cyclotron damping. The Landau and normal 
cyclotron resonant energies for the WWs range from ^1-20 eV and ^25-2100 eV, respectively, 
which supports our intrepretation of the relationship between 8B and T eh ff ec . There is no noticeable 
trend between 8B and T eh /T ec for the LHWs. However, the LHWs were often observed near sharp 
changes in the background magnetic field, suggesting currents may provide some free energy for 
these modes. 

Previous simulation studies [e.g. Saito and Gary 2007 1 examining the effects of whistler mode 
turbulence on suprathermal electrons found that whistler mode waves were capable of scattering 
suprathermal electrons, imposing constraints on temperature anisotropies and strahl width and/or 
density. They imposed wave amplitudes of 5B/B„ ~ 0.07-0.10. The WWs observed herein had 
8BIB„ ~ 0.01-0.16. Therefore, the simulation results support our interpretation of the relationship 
between 8B and |q e |/q and the inverse relationship between 8B and T eh fT ec for the WWs. Further- 
more, these waves are observed primarily when Tu/IV > 1.01 and |q,| ^ 0, which are conditions 
suggested to always be unstable to the whistler heat flux instability \Gary et aT| |1999| . It is dif- 
ficult to say whether the waves are caused by the observed distributions or whether the waves are 
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producing the observed distributions. Independent of the source of either effect, it is clear that the 
separation of the halo from the core is essential to understanding whistler related instabilities in the 
solar wind. 



5.2 Interpretation of Wave Modes 



The polarizations, rest frame frequencies, wave number range, observed proximity to magnetic field 
gradients, and relationship between 5B and the |q c |/q„ are all consistent with the interpretation of the 
LHWs as electromagnetic lower hybrid waves. The positive correlation between 5B and |q P |/q may 
be due to stochastic electron acceleration of the halo electrons, since ©/kii 3> 0)/k ± for these modes 
[e.g. Cairns and McMillan 2005 1. However, below |q e |/q ~ 0.03 we observe no relationship with 



(5B. Perhaps the waves are only affected by the heat flux carrying electrons above some threshold. 
This may suggest that the LHWs, which are often observed near magnetic field gradients, are more 
dependent upon cross-field currents. However, there is no one-to-one correlation between the two 
phenomena. 

The polarization, narrow band frequencies, nearly field-aligned propagation, and bursty wave- 
form features are all consistent with the interpretation of the WWs as electromagnetic whistler mode 
waves. Many of their properties are consistent with the narrow band electromagnetic whistler mode 
emissions called "lion roars," originally observed in the magnetosheath in association with mag- 
netic field depressions [e.g. Smith and Tsurutani 1976 1. Zhang et al. 1 1998 1 observed lion roars 



in the terrestrial magnetosheath not associated with magnetic field depressions and interpreted their 
nearly unidirectional propagation as evidence that they propagated away from some source region. 
\Masood et al. | 2006 1 examined electron distributions looking for relationships between periods sat- 
sifying Equation [T] and observation of lion roars, yet they found no correlation. 

2006| have failed to find a 



We believe that the reason previous studies [e.g. \Masood et al. 



relationship between A,, and lion roars (or whistler mode waves) is due to the use of single bi- 
Maxwellian electron distribution functions. We find that only 5/13 WWs satisfy Equation [T] when 
using the entire distribution, however, that increases to 10/13 WWs when using only the halo. We 
have already shown that all of the WWs satisfy the instability thresholds of Gary et al. | 1994|, which 
requires the separation of the core and the halo. 

Figure [3] shows an example electron distribution that will highlight why we argue for the sepa- 
ration of the core and halo. The values for A,, ,. ;, are shown to the right of the figure. Using these 
values, one can see that only the halo component would satisfy Equation [T] This example clearly 
illustrates why using moments of the entire distribution will not yield values that would predict an 
anisotropy instability. This reiterates our emphasis on the treatment of the electron distribution using 
multiple components when looking for electron-related instability thresholds in the solar wind. 

Whether the WWs are lion roars is beyond the scope of this study. We are interested in how they 
interact with the particles. We do not observe any consistent relationship between the observation 
of WWs and changes in B„ or n e . The correlation between 5B and the |q<.|/q and the fact that all 
the WWs were observed simultaneously with particle distributions unstable to the whistler heat flux 
instability suggest that the free energy for the WWs are the heat flux carrying electrons. Recall that 
10/13 WWs were observed with electron distributions whose halo satisfied Equation [T] In the pre- 
vious section, we interpreted the inverse relationship between 5B and T eh /T ec as increased cyclotron 
damping. Thus, we think the WWs may be causing T ±h /T^, > 1.01 and being driven by a whistler 
heat flux instability. 

The properties of the MIXED modes, at first glance, suggest that these are nothing more than 
the simultaneous observation of LHWs and WWs in the same TDSS event. This can be seen when 
comparing LHWs and MIXED modes satisfying f /o „ in Tableland Figure [T] The MIXED modes 
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satisfying f Mgh show nearly identical propagation characteristics to the WWs, as shown in Section 
4. 1 However, only half of the MIXED modes were observed with electron distributions satisfying 



the instability thresholds of Gary et al. | 1994) and no relationship was observed between SB and the 
|q<.|/q or T eh fT ec . Thus, the MIXED modes do not appear to be interacting with the halo electrons 
in the same way as the WWs and LHWs. 

Since one can see in Figure [T] that all the waves contain some fluctuations at f sc < f a and the 
wavelet transforms show power extending between the low and high frequency components, we 
examined effects arising from a coupling of the two modes. We calculated estimates of threshold 
wave amplitudes for the modulational instability [e.g. Shapiro et al. |1993| and then converted the 
electric field amplitudes to magnetic field amplitudes [e.g. \Huba et al. | ~978| . We find that 10/11 of 
the MIXED modes satisfy the threshold criteria for the modulational instability. However, we also 
find that 18/22 LHWs satisfy the same criteria. Therefore, it is not clear what drives the MIXED 
modes or the LHWs. 



6 Discussion and Conclusions 

We present the properties of electromagnetic lower hybrid and whistler mode waves downstream of 
supercritical IP shocks. The work presented herein can be summarized by the following points: 



1. The lower hybrid wave properties are consistent with theory ^Vlarsch a nd Chan g\\\9%5\\Wu et a 



1983 1 and previous observations at collisionless shocks | Walker et al. 2008 Zhang and Matsumoto 
1998 1. These modes are often observed near magnetic field gradients and they show a very weak 
positive correlation between 8B and the \q e |/q„. However, it is unclear what the source of free en- 
ergy is for these modes; 

2. The whistler mode wave properties are consistent with theory [e.g . \Kennel a nd Petscheck ]\l966\ 
and previous observations \Coroniti et al.\ [T982 ; Len gyel-Frey et al\ \\996\\Zhang a nd Matsumoto 



1998 1. We observe a weak positive correlation between 8B and |q f |/q D and an inverse relationship 



between 8B and T eh /T ec . All of these modes are observed with electron distributions that satisfy the 
whistler heat flux and whistler anisotropy instability thresholds \Gary et arj|1994||1999) ; 

3. The events showing mixed wave features may be due to a modulational instability or they may be 
coincidental observations (i.e. merely simultaneous temporal observations). They show no relation- 
ship between SB and any of the quantities examined for the other two modes. Thus, the source of 
free energy for these modes is unclear; 

4. The separation of the core and halo is absolutely necessary for instability analysis in the solar 



wind. This was clearly demonstrated with our example distribution in Section 5.2 , and 

5. The waves presented herein have amplitudes as large as 8BIB„ ~ 0.3, rest frame frequencies from 

< CO,,, to ~ 0.9D.„, and kp„ ~ 0.01-5.0. 

The observations suggest that the whistler mode waves are receiving free energy from the heat 
flux carrying electrons and producing Ty/Tp > 1.01. The lower hybrid waves may be interacting 
with the heat flux carrying electrons but only for sufficiently large values of |q e |/q . Since lower 
hybrid waves have co/k| ~ W Te 3> (Ofk ± , it may be that smaller heat fluxes do not have noticeable 
effects on the wave amplitudes. At larger values, perhaps these electrons add to the wave amplitude 
or perhaps the larger amplitude waves add to the heat flux. Their proximity to magnetic field gra- 
dients would be consistent with cross-field current instabilities or gradient-drift instabilities, but not 
all are observed near magnetic field gradients. We also observe waves with high and low frequency 
fluctuations independently consistent with whistler and lower hybrid waves, respectively. These 
events have no dependence of 8B on any of the parameters examined in our study. Therefore, the 
only conclusive argument we can make is that the whistler mode waves are driven by the heat flux 
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carrying electrons and cause perpendicular heating of the halo through cyclotron interactions. 

Though we cannot determine whether the waves cause the observed distributions or vice versa, 
we argue that the wave amplitudes clearly interact with the halo electrons. These results have impli- 
cations for theories of the evolution of particle distributions from the Sun to the Earth. For instance, 
the observation of whistler mode waves with unstable electron distributions supports theories sug- 
gesting whistler mode waves limit the halo electron temperature anisotropies and solar wind strahl. 
Note that the observed whistler amplitudes can exceed the estimates used in the simulations [e.g. 
\Saito and Gary||2007| that have led to these theories. These results suggest that the evolution of the 
halo/strahl from the Sun to the Earth could be explained by wave-particle interactions with whistler 
mode waves. However, there is another factor which has not been mentioned. The observed waves 
and distributions all occur downstream of strong IP shocks. This leads to two possibilities about the 
evolution of the electron distributions: (1) the modifications are due to transient events like IP shocks 
or (2) the wave-particle scattering suggested above is a global process but is enhanced downstream 
of IP shocks. Further investigation during quiet solar wind periods is necessary to resolve this issue. 
Future work will compare periods of fast and slow solar wind streams. 
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A Doppler Shift Calculations 



In this section, we define the four methods used to estimate ranges for the wave numbers, |k|, of 
the observed wave modes. Since we are limited to single point measurements, we must make some 
assumptions about the fundamental properties of the observed modes. Below, we will introduce the 
theory and describe how we estimate values for |k|. 

In the solar wind, i2„ <C C0 pe (where CO ps is the plasma frequency and i2„ is the cyclotron fre- 
quency of species s) which limits the number of allowable normal modes of propagation with £l ci 
< CO < £l ce . The three most commonly cited wave modes in this frequency range in the solar wind 
are the magnetosonic \Krauss-Varban and Omidi 1991| , lower hybrid \ Marsch and Chang 1983 1, 



and whistler mode waves | Coroniti et al. 1982| . At low frequencies (f c , < f < f /A ), all three modes 



can lie on the same index of refraction surface [e.g. \Wu et al. 1983). Therefore, let us start with 
the whistler mode. The cold plasma index of refraction \Gurnett and Bhattacharjee , 2005] for an 



obliquely propagating whistler mode wave and CO 2 <C C0 pe 2 is given by: 



2„2 



k z C 

~co~ 2 



2 

CO„e 



co (Qc cos 9 kB - co) 



(2) 



where CO is the rest frame angular frequency and kB is the angle between the wave vector and the 
ambient magnetic field. Solving Equation|2]for |k| will give us our first wave number estimates. 

In the solar wind, the plasma flows with a velocity V„ r relative to the spacecraft frame of ref- 
erence (ignoring spacecraft motion), which causes any observed frequency to be Doppler shifted 
by: 



CO* 



co + k-V,,, 



(3) 
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where co sc is the spacecraft frame or observed frequency. We can solve Equation|2]for co and substi- 
tute into Equation [3] After some algebra, we find: 

= V P + (cos 9 kB - & sc ) P+Vk- & sc (4) 

where k = kcl(O pe , (b sc = (O sc IQ. ce , V = V svl ,cos 9 kv /W Ae , and V Ac = B„/(/i n e m f ) 1//2 . We already know 
9 kB and 9 kv from our MVA results, thus we can solve Equation [4] for |k| giving us our second wave 
number estimates. 

At very low frequencies (co <C £2 ce ) and long wavelengths {kc/co, x <C 1), this new equation 
reduces to the result found by \Coroniti et q/.||1982[ given by: 



k = ± 



V.„,, I cos 9 kv I 
2V Ae I cos 9 kB 



1 ( V sw cos9 kv \ 2 co sc 



(5) 



\2V, e cos 9 kB J £2«. I cos 9 k 



Note that \Coroniti et a/.| ||T982] used the absolute values to indicate that frequencies satisfying co > 
imply (Osc > 0. The physical significance of the ± sign in Equation|5] assuming |cos 9 kv \ = 1, is the 
following: + sign corresponds to a O red-shifted to the observed co sc ; and the — sign corresponds 
to a co blue-shifted to the observed co sc . Results from Equation [5] give us our third wave number 
estimates. 

Alternatively, the cold plasma dispersion relation for an electromagnetic lower hybrid wave is 
given by: 

co\ 2 1 



, M t cos 2 9 kB 

1 - 



1 + co p<; 2 /k 2 c 2 



(6) 



v fl ft y 1 + co„ 2 /k 2 c 2 

where £l lh is given by CO pi /(l + (co p JQ. ce ) 2 ) 1 ^ 2 . Doppler shifting Equation|6]and solving for |k| gives 
us our fourth set of wave number estimates. 

Using the above equations and our observations, we will now present an example calculation 
using the waveform in Figure [T]C. We used two frequency filters on this wave at f sc ~ 5-30 Hz and 
~ 120-200 Hz. Thus, we separated Table[3]into two sections by the different frequency filter ranges 
defined above. The first part of each section of the table shows the background plasma parameters 
(B„, n e , and V w ), filter range, and wave normal angles for this event. Using the range of measured 
parameters shown in the first line of Table[3] we can estimate p„ ~ 1.6 km, cl(O pe ~ 1.3 km, and f !h 
~ 10.3 Hz. Notice that c/co pe ~ p«, thus kc/o ; , f will be ~ kp« for this case. Once we determine |k|, 
we can invert Equation|2]to solve for co using the solutions from Equations [2] [3] and|5]to determine 
the wave phase speed, (o/k. The |k|-values determined from the Doppler shifted Equation [6] were 
put back into Equation [6] to find the corresponding 0)/k. These values compose the last column in 
Table|3] 
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Table 1: Shock Parameters and Mach Number Ratios 



Date 


V shn (km/s) 




M f 


N,. 2 /N ;1 


M f /Mcr 


1998-08-26 
1998-09-24 
2000-02-11 
2000-04-06 


687 ± 27 
772 ± 96 
641 ± 13 
647 ± 98 


82° ± 3° 
82° ± 2° 
87° ± 2° 
70° ± 5° 


4.7 ± 0.2 
2.9 ±0.1 
3.3 ±0.1 
4.0 ± 0.6 


2.9 ± 0.3 

2.2 ± 0.4 

3.3 ±0.5 
3.8 ± 1.3 


2.6 ± 0.4 
1.3 ±0.2 

1.6 ±0.2 

1.7 ±0.3 



Table 2: Summary of Wave Observations 



Wave Type 


#ofTDSS 


# of wave 


<5B/B 




Events 


vectors 


(filtered) 


LHWs 


23 


118 


-0.054 ± 0.003 


WWs 


13 


138 


-0.039 ± 0.002 


MIXED f low 


11 


43 


-0.063 ±0.011 


MIXED f high 


11 


162 


-0.015 ± 0.001 


Total 


47 


461 





16 



Table 3: Plasma parameters and Wave number estimates for Figure [T]C 



Range of Parameters used for Doppler Shift Calculations 


B„ (nT) 


n f (cm 3 ) 


V„ r (km/s) 


f. K (Hz) 


OkB (deg) 


kv (deg) 


-15.7 


-16.0 


-586 


5-30 


56°-57° 


29°-32° 


Wave Numbers from Doppler Shift Calculations 


Method/Equation 


kc/ov 


k(km- 1 ) 


k||P„ 




(o/k (km/s) 


Equa 
Equa 
Equation 
Equation 
Doppler shif 


ti 
ti 

5 
5 
e< 


on 
on 

(- 
(+ 
IE 


2 
4 

s 

5 

q 


ign) 
ign) 

uationpl 


0.1-0.4 
0.1-0.4 
0.06-0.4 
0.3-0.5 
0.1-0.4 


0.2-0.3 
0.1-0.3 
0.05-0.3 
0.2-0.4 
0.1-0.3 


0.1-0.3 
0.09-0.3 
0.04-0.2 
0.2-0.3 
0.09-0.3 


0.2-0.4 
0.1-0.4 
0.07-0.3 
0.3-0.5 
0.1-0.4 


280-660 
270-680 
130-630 
570-820 
280-680 


Range of Parameters used for Doppler Shift Calculations 


B„ (nT) 


n,, (cm 3 ) 


V„ (km/s) 


f sc (Hz) 


6 k B (deg) 


9 kv (deg) 


-15.7 


-16.0 


-586 


120-200 


15°-25° 


82°-89° 


Wave Numbers from Doppler Shift Calculations 


Method/Equation 


kc/cOpt, 


k (km" 1 ) 


k||P„ 




co/k (km/s) 
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Table 4: Wave parameters estimates for each wave mode 



Equation/Wave Type 


kc/(» pe 


kp„ 


k||P„ 


k ± p« 


co/Q. 


Doppler shifted Equation^ 
for LHWs (H = £%) 


0.02 - 5.0 
0.50 ± 0.08 


0.01-5.0 
0.42 ± 0.07 


0.0001 - 3.8 
0.17 ±0.04 


0.003 - 4.6 
0.37 ± 0.06 


0.09 - 5.0 
1.94 ±0.17 


Equation!^] 

for WWs (H = i2«) 


0.2 - 1.0 
0.48 ± 0.02 


0.2 - 0.9 
0.44 ± 0.02 


0.2-0.8 
0.40 ± 0.02 


0.003 - 1.7 
0.15 ±0.01 


0.03 - 0.9 
0.18 ±0.01 


Doppler shifted Equationpl 
for MIXED f /o „ (D. = (O l !^) 


0.03 - 4.9 
0.48 ±0.11 


0.01-5.0 
0.39 ±0.11 


0.0003 - 4.1 
0.15 ±0.06 


0.001 - 4.3 
0.35 ± 0.09 


0.07 - 5.0 
2.01 ± 0.25 


EquationHfor MIXED 

f high (ti = H„) 


0.1 -4.9 
0.92 ± 0.06 


0.07 - 4.9 
0.71 ± 0.05 


0.003 - 4.5 
0.44 ± 0.04 


0.005 - 5.0 
0.50 ± 0.05 


0.003 - 0.96 
0.23 ± 0.01 
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Figure 1: Three examples representative of the three wave types observed in this study including 
the LHW (panel A), WW (panel B), and MIXED (panel C) modes. Panel D shows the right-handed 
coordinate system used to rotate the waveforms. The left-hand column of each wave event panel 
(A-C) shows the FACs of the magnetic field and the right-hand corresponding column the wavelet 
spectrogram. The relative amplitudes of the waveforms are illustrated by the vertical black arrows. 
In the wavelet panels two lines mark the lower-hybrid resonance (black line) and the electron cy- 
clotron frequency (white line). The wavelet spectrograms are uniformly scaled for each individual 
waveform in nT 2 /Hz. On each wavelet spectrogram are a bowl-like line and contours marking the 
cone of influence and 95% confidence levels, respectively | Torrence and Compo, 1998b |. 
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Figure 2: Filtered wave amplitudes vs. T e ,,fT ec and |q P |/q for the WW (panels A and B) and LHW 
(panels C and D) modes. The red squares indicate the mean with the error bars showing the range 
of values for each TDSS event. The heat flux magnitude is normalized by q , which is the free 
streaming heat flux saturation level [e.g. Salem et al. 2003 1 (shown in panel B). The green dashed 
lines are used to show qualitatively the trend for each set of points. 
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Figure 3: An example of an unstable electron distribution observed near a whistler mode wave 
downstream of the 1998-08-26 IP shock ramp. The contour plots show contours of constant phase 
space density in the plane containing the ambient magnetic field (horizontal axis of contours) and 
solar wind velocity. The shock normal direction (red arrow), solar wind velocity (black arrow), and 
sun direction (blue arrow) are projected onto the contour. The bottom plot shows the parallel (solid 
red line) and perpendicular (dashed blue line) cuts of the distribution with associated one-count level 
(solid green line). The vertical black dashed lines show the velocity corresponding to our estimate 
of the spacecraft potential. To the right are values for the right-hand side of Equation[T]for the entire 
(e), core (c), and halo (h) component of the distribution. 
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